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Summary 

The maximum measured electric fields in thunderclouds are an order of magnitude less 
than the fields required for electric breakdown of the air. one explanation for lightning 
initiation in these low fields is that electric breakdown first occurs at the surfaces of 
raindrops where the ambient field is enhanced very locally due to the drop geometry . lab- 
oratory experiments [crabb & latham, 1974] indicate that colliding raindrops which coalesce 
to form elongated water filaments can produce positive corona in ambient fields close to 

THOSE MEASURED IN THUNDERCLOUDS. We CALCULATE THE E-FIELD DISTRIBUTION AROUND A SIMULATED 
COALESCED DROP PAIR AND USE A NUMERICAL MODEL TO STUDY THE POSITIVE CORONA MECHANISMS IN 
DETAIL. OUR RESULTS GIVE GOOD AGREEMENT WITH THE LABORATORY OBSERVATIONS. At THE ALTITUDES 
(AND THUS LOW PRESSURES) AT WHICH LIGHTNING INITIATION IS OBSERVED, OUR RESULTS SHOW THAT 
POSITIVE CORONA CAN OCCUR AT OBSERVED IN-CLOUD E-FIELDS. 



1. Introduction 

Lightning initiation in thunderclouds is poorly understood. An order of magnitude 
discrepancy exists between the maximum measured electric fields (E-fields) in clouds 
(E m ax) and the E-ficlds required for dielectric breakdown of air. E max ~ 100 - 400 kV/m 
Marshall ct al, 1995| , Winn ct al, 1974 . Ebreakdown ~ 2700 kV/m at surface pressure; at 
the lower pressures (~ 500 mb) at which most lightning is observed to initiate, Ebreakdown 



is reduced to ~ 1400 kV/m - still far greater than E max . 

One of several explanations put forward to explain this discrepancy is the enhance- 
ment of the electric field near the surfaces of hydrometeors (water or ice particles) in 
clouds. A set of laboratory experiments by Crabb & Latham (hereafter CL) showed that 
colliding raindrops may provide the starting point for lightning initiation. CL obtained 
very promising results in a set of experiments in which they measured the E-fields re- 
quired to initiate a discharge from the surface of filamentary, temporarily coalesced drops 
created when two water drops collided. They observed pulsed, intermittent discharges in 
a localized region near surface of the drop and found that the E-fields required lay within 
the range of observed thunderstorm E-fields. 



We extended models of the positive discharge process developed by [Dawson & Winn, 1965 
Gallimberti, 1979, Bondiou & Gallimberti, 1994, Abdel-Salam et al, 1976[ in order to study 



the discharge processes occurring from hydrometeors. CL's laboratory conditions were 
used to initialize the model and their results were used to validate it. With our model 
we were able to vary both the microphysical and environmental conditions and inves- 
tigate a range of conditions applicable to those found in thunderclouds. In particular, 
we investigated continuous discharges from the drop surface and we studied the pressure 
dependence of discharge initiation E-fields. 

We begin with a brief description of Crabb & Latham's experimental procedures and 
results, followed by a discussion of the discharge model used to make the calculations. 
Finally we discuss our model results - showing the E-field required to initiate various 
discharge types as a function of the coalesced drop properties and air pressure. 

* Corresponding author: Geophysics Program, University of Washington, Box 351650, Seattle, 98195- 
1650, USA email:vicki@geophys. washington.edu 
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2. Laboratory Experiments 



(a) (b) 




Figure 1. (a) Schematic of Crabb & Latham's experimental setup in which two drops (R=2.7mm and 
r=0.65mm) collided in the presence of an applied electric field, (b) Photograph of the coalesced drop 

that formed after the collision. 

Figure]]] (a) shows a schematic of the CL experiment . Their chamber, held at surface 
pressure, had a positive, high voltage upper plate and a grounded lower plate separated 
by 50mm. Voltages of up to 30 kV could be applied - corresponding to a maximum 
uniform E-field of 600 kV/m within the chamber. Large water drops (R=2.7mm) were 
dropped into the chamber and collided with small drops (r=0.65mm) which were blown 
upwards, simulating drops moving in updrafts in thunderclouds. A variety of coalesced 
drop shapes were observed, depending on the nature of the collision. CL described three 
basic collision modes: head-on, glancing and intermediate. Glancing and intermediate 
collisions produced a coalesced drop with a long filament extending from the large drop 
- see Fig [j] (b). Head-on collisions resulted in a flattening of the large drop and did not 
produce these long filaments. The drops remained in the coalesced state for ~ 1 ms. 

In CL's setup, a negative charge was induced on the upper surface of the drop while 
the lower end had a positive induced charge. In the thundercloud setting these drops 
would be located above the negative charge center of the cloud. CL recorded the size and 
shape of the coalesced drops as well as the applied E- fields required to initiate discharges 
for a large number of coalesced drops. They observed discharges from both ends of the 
drop but focussed on the positive pulses occurring at the lower surface of the drop. This 
surface was observed to remain intact. In contrast surface disruption was observed at the 
upper, negative surface of the drop. 

CL observed that positive burst pulses occured for values of E between 250 and 500 
kV/m, depending on the length of the coalesced drop. 

3. Model Description 

We describe three basic discharge processes that can occur at the surfaces of drops 
in the presence of strong E- fields. 
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The first process, surface disruption discharge, occurs when the electrostatic repul- 
sive force on a drop in a strong E-field exceeds the surface tension. This results in breakup 



of the drop surface, and an associated discharge. Dawson, 1969 observed the surface E- 
field, Edisruption, required to initiate this form of discharge as a function of drop size. 
Edisruption is independent of pressure. 

The remaining processes are referred to as 'pure' corona processes because the dis- 
charge initiates without the occurrence of drop surface disruption. These processes are: 

• burst pulse discharges, which are intermittent, and 

• continuous streamers, which are capable of propagating continuously. 

We discuss each in detail below. All of these processes result in deposition of charge 
on the drop: either positive or negative charge depending on the sign of E exter nai- We 
focussed only on positive corona: it is simpler to model, has a lower initiation threshold 
and was studied more extensively by CL than negative corona. 

(a) Positive Pure Corona model 




Figure 2. Schematic of positive discharge formation near the surface of a drop. Free electrons are 
accelerated by the E-field and undergo collisions with air molecules. The ionization of molecules within 
the ionization zone leads to an exponential growth of electrons (avalanche) and the formation of a 

spherical streamer head. 



Following pawson fc Winn, 1964 |Gallimbcrti, 1979| , |Bondiou fc Gallimbcrti, 1994 
we model the positive discharge as a series of electron avalanches. 

Consider the electric field near the surface of a drop which is situated in an external 
electric field, E externa i. (See Fig |[) The total electric field is a function of the distance, 
r, from the drop surface. Initially, the total electric field at r is 



E(r)=EJr)=E, 



external 



+ Edrop(r) 



(1) 
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where Edrop is the contribution due to charge induced on the drop. E g {r) is called the 
geometric field. 

In the presence of E free electrons are accelerated and undergo collisions with air 
molecules. At some radial distance from the drop E is such that: 



a(E/p) = r)(E/p) 



(2) 



where a [m _1 ] and rj [m _1 ] are the ionization and attachment coefficient for electrons 
in air, respectively and p is the total air pressure. The surface defined by eqn (^) is the 
ionization zone boundary - inside this boundary a > i] and there is a net growth of free 
electrons. At surface pressure the ionization zone boundary is the surface along which 
2700 kV/m . Figure |3| shows a and rj as functions of E/p jGeballe, 1953 Loeb, 1965 



E 



Badaloni, 1972, |Ibrahim fc Singer, 19821 



Following [ pawson fc Winn, 1964 [Griffiths fc Phelps, 1976| , pallimbcrti, 1979 



we 

replace the three-dimensional problem by a one-dimensional one in which all avalanches 
occur along the z-axis. The point Zi marks the intersection of the ionization boundary 
with the z-axis. When a free electron, starting at Zj, is accelerated by E towards the 
drop, the number of electrons grows exponentially as z decreases. This is referred to as 
the primary electron avalanche. Due to the exponential nature of the growth, most of 
the ionizing collisions occur near the surface of the drop. The free electrons are then 
absorbed by the drop, leaving behind a concentration of positive ions, modeled as a 
sphere Dawson fc Winn, 1964 pallimbcrti, 1979 1 , and referred to as the streamer head. 

The number of positive ions formed by the primary avalanche traveling from the 
ionization zone boundary, Zi, to the drop surface, z 0l is given by: 



(3) 



Ni = exp[ / {ct(z)) — j](z))dz] 

J Zi 

The radius of the streamer head is approximately: 

o 

R s = {6 [(D(z)/v(z)dz)} 1 / 2 



(4) 



where D and v are the electron diffusivity and drift velocity, respectively. D and v are 



E{z) = E g {z) + 



47re (z - Rs) 2 



(5) 



where the second term, E c , is the E- field due to the spherical charge concentration of the 
streamer head. 

In addition to ionization, collisions between the free electrons and air molecules also 
result in the excitation of the molecules, which then emit photons on decay. A certain 
fraction of these photons in turn have sufficient energy to ionize molecules that they 
encounter, creating photoelectrons. These photoelectrons then start a series of secondary 
avalanches which converge on the drop from all directions. 

The number of photoelectrons created per m at a radial distance, from the drop 
surface is given by: 



functions of the ratio (E/p) and thus depend on z [Healey & Reed, 1941, Ibrahim & Singer, 1982] 
R s « 30/xro at surface pressure for most of the calculations reported here. 
The total electric field at z is now given by: 



P(i) = AJVi • expl-^1] -h-G 



(6) 
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(a) Ionization coefficient 



(b) Attachment coefficient 
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Figure 3. (a) Ratio of ionization coefficient to pressure, a/p, for electrons in air [ Badaloni. 1972 



eballe, 
IT 



Loch, 1965 1 fb ) Ratio of attachment coefficient to pressure, rj/p, for electrons in air [Badaloni, iy/2, 
"T953I (c) = fi ■ f-2 -8 where /i is the number of photons created per ionizmg collision, fi 
the number of photoions crea ted per photon per meter and 9 is a solid angle = 2-7T in our 
calculations [Penney & Hummert, 197CL (dj Batin nf phntnn ahsprptinn coefficient to pressure, p./p, in 

air [Penney & Hummert, 197C]. 
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where /i is the number of photons created per ionizing collision 
H [m _1 ] is the photon absorption coefficient in air 
/a [m _1 ] is the number of photoions created per photon per meter 
G is a geometric factor to account for the fact that some 
photons are absorbed by the drop. 

Both n and /i • j% are functions of I ■ p, the product of the distance from the photon 
source (the collisions) and air pressure | Penney fc Hummcrt, 1970 - see Fig || 



Then the total number of ions created in the secondary avalanches is given by: 

N 2 = I ° P(l) ■ exp[ f \a-rj) dz] dl (7) 



where z Q indicates the position of the primary streamer head surface. 

(i) Initiation Conditions 

A burst pulse discharge is initiated if the number of photoelectrons created along 
the ionization zone boundary during the growth of the primary avalanche is equivalent 
to the number of photoelectrons that started the primary avalanche (commonly taken as 



1) [Abdel-Salam et al, 1976 



We consider photoelectron production in a region of depth (l//x) along the ionization 
zone boundary and write the above condition as follows: 

S2> 1 (8) 

This type of discharge is intermittent because the number of positive ions, N±, in the 
primary streamer head is too small to attract the following avalanches to its surface. 
Instead, the successor avalanches are directed to the drop - allowing the discharge to 
"spread" over the drop surface. 

The fulfillment of eqn (||) is strongly influenced by the relationship between the 
mean free photon path (l//x) in air and the location of Zi. For drops with small radii, Zi 
is closer to the drop surface than for those with larger radii. Thus for small drops the 
number of photons likely to reach the ionization zone boundary is increased. The chances 
of a photoelectron being produced then increases and it is easier to initiate a burst pulse 
discharge under these conditions. 

A more stringent initiation condition exists for continuous streamers. In this case 
the number of positive ions in the primary streamer head must be large enough to attract 
the secondary avalanches to the streamer head surface. This is achieved when the radial 



E-field around the streamer head, E c — AlTC ^l R p ~ E g [ Abdel-Salam et al, 197f ]. In 
addition: 

(a) N2, the number of positive ions in the streamer head that results from the secondary 
avalanches, must equal Nx, the number of positive ions created by the primary avalanche, 
and 

(b) the radius of the secondary streamer head must equal R s , the radius of the primary 
streamer head. 

These conditions ensure that the initial streamer head charge density is reproduced in 
the second streamer head. Continued reproduction of the streamer head in subsequent 
steps results in propagation of the positive streamer away from the drop surface. For all 
the geometric conditions considered in this paper, the initiation of continuous streamers 
requires a larger external E-field than for burst pulse discharge initiation. 



CORONA EMISSION FROM HYDROMETEORS 



7 



The minimum value of E ex temal necessary to initiate a discharge at pressure p 
is referred to as Ei n iti a tion(p) and depends on the type of discharge (burst pulse or 
continuous streamer). 

(6) Model Procedure 
Our results were obtained using the following procedure: 

• Define the drop shape and permittivity, e. The idealized shape used is shown in Fig||[ 
Set the air pressure, p. 

• Apply Eextemal to the drop. 

• Calculate the E-field distribution around the drop using a finite element method based 
solving routine [Quickfield]. 

• Compare the E-field at the drop's negative surface to the known surface disruption 
field threshold, E 'disruption [ Dawson, 1969 ], 

If E(surface) > E 'disruption then add varying amounts of positive charge, Qdropy to the 
drop. 



Recalculate the E-field distribution [ Quickfield |. 
Find the ionization zone boundary, Zj. 
Compute Ni and R\ from eqns (||) and (||) respectively. 
Compute P(l) at Zi from eqn (g). 

For = 1, E ex ternat = E ini tiation{p) for burst pulse discharges. 



Compute N 2 and R 2 from eqns (0) and (|J) respectively. 



EE r 
ers. 



E„,N 2 = iVi and R 2 = Ri, then E, 



external 



Ej. 



t (p) for continuous stream- 



■external 




negative discharge 



idealized model 
drop 



positive discharge 



Figure 4. Idealized, model drop used in this study with lines of equipotential shown. Th e field distri- 
bution around the drop is calculated using a finite element method [Quickfield |. 
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4. Results 



(a) Surface disruption 
For E ext emai > 200 kV/m the calculated E-field at the surface of the upper, negative 
end of the drop was > 8500 kV/m, th e value of Edi srU ption at p=1000mb for a water drop 
of radius r =0.65mm |Dawson, 1969|] . This is consistent with CL's observations that the 
upper surface disrupted at these E-field strengths. The resulting negative discharge then 
deposited positive charge on the drop. 



(6J Einitiation VS Qdrop 

Fig |^ shows the Einitiation values for positive burst pulse discharges from the lower 
positive end of the drop as a function of the charge, Qdrop, deposited on the drop by the 
negative discharge from the upper end. The drop length is held fixed at L=20mm. 




Figure 5. Ei n i t i a tion for positive burst pulse discharges from the lower positive end of the drop as a 
function of Qdrop? the charge deposited on the drop by the negative corona from the upper end. The 

drop length is held fixed at L=20mm. 

Einit iation decreases rapidly once Qdrop exceeds 10 -10 C. The Rayleigh stability 



criterion | Rayleigh, 1882 , Taylor, 1964 gives Qrl, the maximum charge that a sphere of 
liquid can hold before the electrostatic repulsive force overcomes the surface tension. In 
SI units it is given by: 

Qrl = 64 ir 2 e a -r 3 a (9) 

where r is the sphere radius and a is the surface tension. 

For our drop dimensions Qrl ~ 4 x 10 -9 C. Since CL did not observe disruption 
of the lower surface of the drop, we limited our calculations to Qdrop < Qrl- For larger 
allowed values of Qdrop, close to the Rayleigh limit Qrl, the values of Equation become 
comparable to CL's experimental values and to those observed in thunderclouds. 
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In addition to the burst pulse discharges we also calculated the fields required to 
initiate continuous streamers. For Qdrop just below the Rayleigh limit, Equation ~ 400 
kV/m for continuous streamers, approximately 50% greater than that required for burst 
pulse discharges. 
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pi gure 6. ^initiation for burst pulse discharge as a function of the drop length L, for fixed charge 
density. Triangles: calculated values of Equation for burst pulse discharges. Circles: Crabb & Latham's 

measured values. 

(c) E init i ation vs drop length, L 

We now held the charge density, p, on the drop fixed at 0.035 C/m 3 and varied the 
drop length, L. The circles in Fig ^ represent CL's measured values. We found that our 
modeled values of E initiation for the burst pulse discharges (A) decreased with increasing 
L, consistent with the trend that CL observed. 

The agreement between the calculated results and observation is promising and offers 
validation of our model processes. The scatter in CL's results is most likely due to cither 
the differences in the shape of the lower end of the coalesced drops or the amount of 
charge that is deposited by the negative discharge. We found higher Einitiation values for 
a coalesced drop with a spherically shaped lower end, while lower Einitiation values were 
recorded for more pointed lower ends. Our idealized shape with p=0.035 C/m 3 , however, 
provided good agreement with CL's average values for Einitiation- 

The same calculations were carried out for continuous streamers and the results are 
shown in Fig |?]. Einitiation decreased with L in much the same way as for burst pulses. 
The Einitiation values for the continuous streamers were, however, ~ 50% larger than 
those required for burst pulses. 

Figure || shows the two competing processes that determine the dependence of 
Einitiation on L for continuous streamers. On the one hand, for a given ambient E- 
ficld, the surface field at the tip of the filament increases with increasing L, which lowers 
Einitiation- In opposition to this, as L increases E g (z) decreases more rapidly with z, 
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L[mm] 10 400 pressure [mbar] 

Figure 7. Ei„iti a tion for continuous discharges as a function of both drop length, L [mm] , and pressure, 

p [mb] 



the distance from the surface. This reduces the size of the ionization zone and thus in- 
creases Einitiation- Fig (7| shows that the former process dominates; i.e. that the increased 
average field within the ionization zone compensates for the electron's shortened path - 
leading to a lowering of Emulation as the filament length is increased. As Fig indicates, 
dEinitiation / dL decreases as L increases and the effect of increased length becomes less 
significant for L>20mm. 




Figure 8. E-field as a function of distance, z, from the surface of the drop. E s ; and E ss are the surface 
fields for long and short drops, respectively. The ionization zone boundaries for long and short drops are 

indicated by zu and Zi s , respectively. 
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(d) The pressure effect 

All CL's measurements were made at surface pressure (1000 mb). It is, however, of 
interest to know what the Emulation values for continuous streamers would be at the 
lower pressures found in the regions where lightning initiates. We therefore calculated 
Einitiation for continuous streamers for a range of pressures. 

The variation of Equation for continuous streamers with both pressure and drop 
size is shown in Fig ^. The dark region in the lower left corner indicates the region in 
which initiation is most favorable - large L and low pressure. Over the chosen ranges of 
pressure and L, pressure has a greater effect on E initiation than L. 

Fig |t] indicates that E initiation varies linearly with pressure. We consider the de- 
pendence of the various parameters used by the model: a, r/, D and v are functions of 
E/p while the /i and f\ ■ ji are functions oil ■ p. The linear dependence of Equation for 
continuous streamers suggests that the dependence on E/p dominates and that there is 
a unique value of the "reduced" E-field, Y initiation = E initiation / V for a particular E and 
p combination. 



(e) Propagation 

The E-field necessary to sustain stable streamer propagation, E propaga tion{p), was 
measured by [ Griffiths fc Phelps, 1976b as a function of air pressure, p. These stable 
streamers, once initiated, will continue to propagate provided Einitiation > E vrovagat i n- 



Griffiths and Phelps found that E propagat ion 
that E, 



propagation 



(p) oc p 5 (Fig |). At p = 500 mb E. 



400 kV/m for dry air at p = 1000 mb and 



propagation 



150 kV/m for dry air. 
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Figure 9. 



gation as a function of absolute humidity and pressure, p Griffiths & Phelps, 19761 ] 



At p = lOOOmfr Fig shows that Equation > 400 kV/m for all L. These initiated 
streamers will therefore be able to propagate over the entire length of the region in which 
Eextemai remains constant. In thunderclouds this scale is typically hundreds of meters. 



At lower pressures E ini 



>E; 



propagation 



over a large range of L. 
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5. Discussion 

In this paper we have shown that continuous, propagating streamers can be initi- 
ated from water drops at E-fields found in thunderstorms. Provided that Ei n i t i a uon{p) > 
Epropagation (p) , these streamers are capable of propagating over considerable distances. 
This distance is limited by the size of the region in which E exter nai is greater than 

^propagation 

When the electron currents in streamers become large enough, Joule heating pro- 
duces a 'warm' leader; a channel in which thermodynamic equilibrium is destroyed and 
hydrodynamic effects become important. This is commonly referred to as the 'stepped- 
lcadcr' in the cloud-to-ground lightning context. The currents carried by individual 
streamers initiated at the drops are several orders of magnitude too low to produce 



leaders [Bondiou, 1997 1 . These streamers may, however, still eventually lead to leaders if 
they can be combined or multiplied. 

Griffiths & Phelps [1976] considered the role of small scale discharges in thunder- 
clouds, calculating the E-field enhancement due to multiple propagations of positive 
streamers near an electrode. According to their model, a series of three to seven stream- 
ers gave rise to an enhanced E-field of up to ~ 1500 kV/m in a region of several meters 
in linear scale near the electrode. It is possible that several continuous streamers ini- 
tiated from drops in the thundercloud could provide the required field enhancement. 
| Griffiths fc Phelps, 19761 bund that the field was intensified on a time scale of ~ 1 ms, 



which is comparable to the lifetime of the coalesced drops as measured by CL. Further 
investigation is required to determine whether a hydrometeor is capable of initiating 
multiple streamers. 

An alternate mechanism for leader formation would be the combination of several 
streamers in close proximity to one another to form a single, more vigorous streamer with 
sufficient current to transform it to the "warm" leader stage. If we think of drops that 
initiate continuous streamers as "electrodes" then the number of "electrodes" available 
increases with increasing E (see Fig ^|). Thus the likelihood of several streamers initiating 
in close proximity increases and the chance of leader formation is increased. This is also 
in keeping with the observations of large amounts of corona activity in thunderstorms 
without lightning. Only if the "electrode" density is sufficiently high will streamers be 
able to merge and form a leader. These possible mechanisms for leader formation require 
further investigation. 

The streamers observed by both CL and examined in our model were all positive, 
occuring at the lower end of drops. This corresponds to drops located above the neg- 
ative charge center in clouds. Leader formation in this region is observed to lead to 
intra-cloud lightning flashes. Drops located below the negative charge center have neg- 
atively charged lower ends and investigation of this situation will require the modeling 
of negative streamers which are much more complex in nature than positive streamers 



Castellani et al, 1994 . Leader formation in this region (below the negative charge cen- 
ter) will lead to cloud to ground lightning. No attempt has been made in this paper to 
model these negative processes but future attempts should investigate this phenomena. 

Finally, while we have concentrated on liquid hydromctcors, future work should 
incorporate ice particles as possible "electrodes" . This may explain how lightning initia- 
tion can occur at higher altitudes near the upper positive charge center in thunderclouds 
where there is little or no liquid water available. 
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